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Abstract 

The effect of the electrodynamic forces on a charged particle in 
a propagating plane electromagnetic wave is investigated. First it is 
pointed out that for constant fields fulfilling the radiation condition 
there will be an acceleration in the direction of the Poynting vector. 
When oscillating fields are considered the Lorentz force on the particle 
only causes a drift, with constant average velocity, in the direction 
of propagation of the wave, i.e. the direction of the Poynting vector. 
Finally, when the radiative reaction (radiation damping) force is added 
the result is again an acceleration in the direction of wave propagation. 

PACS classification numbers: 03.50.De, 41.60.-m, 41.75.Jv 
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1 Introduction 



Recently Professor Thomas Gold pQ has published a manuscript on the Web 
stating that radiation pressure from the Sun does not exist, or at least that 
it cannot be used for propulsion with solar sails, as has been suggested [2] . I 
do not understand Gold's, mainly thermodynamic, arguments but the phe- 
nomenon of radiation pressure is well established theoretically and exper- 
imentally so the manuscript is somewhat surprising. Usually the force of 
radiation pressure is explained by discussing how photons carry momentum 
which is either absorbed or reflected when impinging on a body. Here we will 
point out and elucidate the less known fact that the electrodynamic forces 
(Lorentz force and radiation damping force) on a charged particle from a 
plane electromagnetic wave accelerates the particle in the direction of prop- 
agation of the wave. This, at least, demonstrates that radiation pressure is 
an immediate consequence of the relativistic equation of motion of a charged 
particle in an electromagnetic field. 

For simplicity of notation we will use Gaussian units. The relativistic 
equation of motion for a charged particle of charge q and mass m, in an 
external electromagnetic field F a b = (E, B) (with the notation of Landau 
and Lifshitz 3j) is given by, 

du a _q b 
m—--F b u, (1) 

if we, for the time being, neglect radiative reaction. Here u a = (1/ c)dx a / dr = 
7(1, v/c) = (u°, u), where x a = (ct, r), and r is proper time, dr = ^dt, and, 
finally, 7 = u° = 1/ \Jl — v 2 /c 2 . Later we will add g a , the radiation damping 
force, on the right hand side. This is a correction needed since an accelerated 
charge radiates and this results in a reaction force. If one introduces the new 
independent variable 

C = SL, (2) 

mc 

one can write Eq. ((!} in the form 

= (u ■ E, 7.E + u x B) (3) 
for the time and space components respectively. 
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2 Motion due to the Lorentz force of constant 
fields 



The general solution of Eq. (j3J) for constant electromagnetic fields F ab = 
(E, B) is known and has been discussed extensively in the literature. In 
particular one can recommend the studies by Salingaros jU |3] . Other illumi- 
nating contributions are by Hyman [Hj and by Munos [Zj. 

First we note that the case of radiation is very special and differs from 
the general case. For constant E and B one can in general make a Lorentz 
transformation to a reference frame in which these two vectors are parallel. 
Then the Poynting vector, 



which represents the flux density of momentum in the electromagnetic field, 
becomes zero. In that frame therefore there can be no radiation pressure 
force. In such general fields charged particles will have a drift velocity equal 
to the velocity of the frame in which the fields are parallel. Radiation fields 
on the other hand are characterized by 



and for such fields there is no such reference frame of zero Poynting vector. 
The general solution of Eq. (jSJ) for constant E and B fulfilling the radiation 
conditions, for a particle starting at rest, is given by Salingaros [H] and is 



The meaning of this result is that, for small ( the velocity will be essentially 
parallel (or anti-parallel) to the electric field as intuition demands, but for 
large ( the velocity will become more and more parallel (never anti-parallel) 
to the Poynting vector S = (c/4:7r)E x B. 

3 Motion due to the Lorentz force of of prop- 
agating wave 

The result (J7j) by Salingaros and its interpretation does not survive when the 
field is allowed to have an oscillatory time dependence, as we will now show. 



E-B = 0, 

\E\ = \B\, 



(5) 
(6) 




(7) 
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We start by rewriting Eq. (|T|). 



du a 



on matrix form, as follows 



dC = F(C)u ' 



(8) 



(9) 



and we will allow F to depend on time via (. Here the components of the 
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(10) 



respectively. We now specialize to the radiation case and chose the x-axis 
in the direction of E so that E = E(()e x . If also chose the the 
direction of the Poynting vector we must have B along the y-axis and of 
the same length as E, so that B = E(()e y . The Poynting vector is then 
S = ^-E 2 (()e z , and the F-matrix is given by 



/ E(C) 

E(0 -E(() 



V E{() 

where we have defined the matrix 



E(()H, 
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'12) 



We have now found that, in the radiation case, the equation of motion, Eq. 
(JOJ), becomes 

^ = [£(C)H]u. (13) 

Thus, differentiating u with respect to ( multiplies u with a matrix. The 
general solution of this equation is given by 



u(C) = exp 



E(r])dri H 



u . 



(14) 



4 



To get an explicit solution we put 



/(C) 



E(r))dr), 



(15) 



and use the series expansion of the exponential. For this we need the powers 
of the matrix H and these are: H° = 1, the four by four unit matrix, H 1 = H, 



/ 1 -1 \ 





V i o o -l / 



(16) 



All higher powers of H are zero matrices. We thus find that 

1 



"(0 



1 + /(C)H + -/ 2 (C)H 2 



u . 



(17) 



We now assume that the initial condition is Uo 
particle is at rest. Explicit calculation then gives 



[1 0), i.e. that the 



u(C) 



/ i + |/ 2 (C) \ 
/(C) 
o 

V |/ 2 (C) / 



If we take E(() = const, we recover Salingaros' solution (JZJ). 

If we instead assume that we have a simple harmonic wave, so that E(Q = 
Ex cos(u>£) and f(() = E x sm(w()/w. We then find that the time average of 
the four velocity (fTHj) becomes 



<u( C )> cS - J ( u(0d C 



V 





hr 2 



(19) 



The acceleration in the direction of the Poynting vector (the z-axis) that we 
found in the case of constant fields has now become a drift with constant 
average speed (originally calculated by McMillan |H]). The speed of this drift 
is larger the smaller the frequency w is, but there is no acceleration and 
thus no average force in the direction of the Poynting vector. For a recent 
discussion of this problem, see McDonald and Shmakov [9 . 
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4 The radiative reaction force from propagat- 
ing plane wave 

So far we have neglected the radiative reaction force and used the equation 
of motion, Eq. To make it more accurate we must add the four force g a 
on the right hand side. Now we investigate the form of g a and how it will 
modify the solutions found in the previous sections. 

The origin of this force is the electromagnetic radiation that an acceler- 
ated charged particle sends out. This radiation carries energy and momentum 
and there is therefor a reaction force on the particle itself. One can show 
that the force, due to dipole radiation, should be / 



V 

3c 3 



v. The four vec- 



tor form of this should give the four force g$ = f^^pr- A four force must 
for kinematic reasons fulfill g a u a = 0, the four scalar product with the four 
velocity must be zero. This is achieved by modifying by subtracting the 
component along the four velocity: g1 = gQ — (gQUb)u a . Since u a u a = 1 this 
gives the desired property. Finally one can insert the expression for 
obtained by differentiating Eq. (JTjl, and also the expression for from the 
same equation. The force g a should, after all, be a perturbation compared 
to the Lorentz force. Thus one obtains 



(20) 




(21) 



Some algebra, and use of the fact that F ab u a Ub = (due to the anti-symmetry 
of F ab ), turns this into 



9" 



2q 3 dF a b 



(F" C F% - F d c F\u iU «) u b . (22) 



3mc 4 dr 3m 2 c 5 
The corresponding three dimensional force / can be obtained by noting that, 
with our conventions, g a = (j/c)(f ■ v/c, f). 

Landau and Lifshitz j3] (§76, Problem 2) give the three dimensional form 
for this force as follows 



4 q 



3 m 2 c 4 



2 q 3 fdE v dB\ 



E x B ■ ( - x Bj x B+ ■ E) E 



(23) 
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Here ^ = + v ■ V and 7 = 1/y 1 — v 2 /c 2 . We will now analyze the effect 
of this force for a charge in an oscillatory plane wave. 

The first term, f , in line ()23j) is then of an oscillatory character. It can, 
in fact, be combined with the Lorentz force, 



v 



f l (t) + f (t) = [qE(t) + q-xB(t) 



2 q 3 fdE v dB\ 

^ dF + c x dT (25) 



3 mc 



v 



qE(t + e t )+q- x B(t + e t ) 

c 



fi(t + e t ), (26) 



and we see that the effect of this term is to make a tiny time translation of 
the Lorentz force with e t = §£37. This essentially is the time needed for 
light to pass across the classical electron radius. This is unlikely to have any 
noticeable physical effects. 

The essential radiation damping force is thus f s . All its terms ()24j) are 
quadratic in the fields so there is hope that this force, / a , produces some net 
work on the particle. In the rest frame (v = 0) we simply get, 



(27) 



and therefore the particle starts to accelerate in the direction of the Poynting 
vector. What happens when it picks up speed? 

Assume that the field is of the type in Eq. (jllj) with oscillatory E = E((). 
Since E = Ee x , B = Ee y , and Ex B = E 2 e z , some algebra then shows 
that our radiation reaction force, quadratic in the fields, is 



2 q 4 



3 m 2 c 4 



E 2 [l 



7 2 d 



:28i 



We see that the effect of this force is always to accelerate the particle in the 
direction, e z , of the Ponting vector. The dissipative, —v, term is one order 
higher in v/c and will thus only become important near relativistic speeds. 
The force goes to zero when v — > ce z , but it does not change direction. 

Surprisingly we find that this force, which usually is called the radiation 
damping force, does not damp anything in this case. Instead its effect is 
the acceleration of a charged particle in the direction of propagation of a 
travelling electromagnetic wave. 

5 Conclusions 



Above we have presented an elegant matrix method that yields an exact 
solution for the motion of a particle in an electromagnetic radiation field with 
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a given time dependence. This generalizes some older results for constant 
fields of Salingaros |UE]. The concise explicit formula of Eq. ()28|) for the 
radiation damping force in a plane wave might be also be a new contribution. 
The main aim of this article has been to illuminate the way in which classical 
electrodynamics leads to a pushing of charges along the direction of of wave 
propagation. 

Just as in solar sailing, this is a very small force, but in space where 
there normally are no dissipative forces, even a small force can have very 
large effect in the long run. Note that there is no need for a force to be 
stronger than gravity to accelerate things outwards from the sun (a common 
misunderstanding). A typical object orbits the sun in a stationary ellipse. 
When a small force is added this ellipse will gradually change, and with time, 
these changes can become very large. 
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